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Abstract 

We show that the evolution equation of the effective potential in 
the auxiliary mass method corresponds to a leading approximation of 
a certain series. This series is derived from an evolution equation of an 
effective action using a derivative expansion. We derived an expression 
of the next-to-leading approximation of the evolution equation, which 
is a simultaneous partial differential equation. 
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1 Introduction 



Finite temperature field theory, which is based only on the statistical mech- 
anism, adequately describes many physical phenomena such as phase transi- 
tions and mass spectrum in a thermal bathjl], |, [J. Its perturbation theory, 
however, breaks down at a high finite temperature. Mass squared becomes 
negative even in a vacuum at the finite temperature. This problem is solved 
by the ring resummation, which adds thermal mass to zero temperature mass 
beforehand^, [5|]. This procedure is still insufficient yet to make the pertur- 
bative expansion be reliable, especially around the critical temperature. For 
example, the perturbation theory indicates that the phase transition of the 
Z 2 -invariant scalar theory is of first order incorrectly ||. For the other exam- 
ple, different properties are indicated by the perturbation theory and lattice 
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and the Standard model [U, 0, p|, [T|, [20 
large Higgs boson mass range. 

These failure of the perturbation theory is caused by bad infrared behavior 
around the critical temperature P5| , . A loop expansion parameter becomes 
XT/m even after the ring resummation due to this infrared effect. Here A 
is a small coupling constant and m is mass at the temperature T. The 
perturbation theory is, therefore, unreliable at a high temperature and small 
mass. This situation arises around the critical temperature of second order 
phase transitions or weakly first order phase transitions. The auxiliary mass 



method controls this infrared behavior by introducing "auxiliary mass" [[27 



53]. We first calculate an effective potential with the large 
auxiliary mass, M ~ T. This effective potential is reliable since the loop 
expansion parameter is small thanks to the auxiliary mass. We next calculate 
the effective potential at true mass from this effective potential through an 
evolution equation. We solve the evolution equation of the effective potential 
with respect to variation of mass squared. 

We used a certain approximation to derive the evolution equation in 



Ref . pq, |29| . Though we got quite good results using the approximation, we 
did not have methods to improve the approximation. In the present paper, 
we show that the previous evolution equation is a leading approximation of a 
certain series. We then derive a next-to-leading evolution equation, which is 
a simultaneous partial differential equation. Though it is difficult to solve the 
evolution equation at an arbitrary temperature due to numerical problems, 
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we can, in principle, improve the approximation systematically. 



2 Evolution Equation 

In this section, we explain our idea and derive the evolution equation for the 
^-invariant scalar theory using the auxiliary mass method. 1 Let's consider 
the following Euclidean Lagrangian density with mass squared m 2 , 

Ce& m 2 ) = -~ (0) 2 - ^ W) 2 - ~m 2 2 - ±{<?f + J m <P + c.t. . (1) 

Here J m {x) is an external source function which sets 

< >= C ( X ) (2) 

and hence depends on m 2 . We set true mass squared negative, m 2 = —fi 2 , 
since we investigate the phase transition of this theory. We assume that the 
coupling constant A is small so that the perturbation theory be reliable at a 
low temperature (T <C m) . 

Our idea is the following. The effective action for the theory, r[<j) c ;m 2 ), 
satisfies the following identity, 

T[0 C ; V) = /J 2 {P^ 1 ) dm' + T[0 C ; M 2 ). (3) 

We set M 2 so large (~ T 2 ) that we can calculate the initial condition, 
T[0 C ; M 2 ), reliably by the perturbation theory. If we can evaluate the deriva- 
tive dT((f) c ) / dm 2 correctly, we can calculate the effective action accurately^. 



2.1 Derivative of the effective action with respect to 
mass squared 

In this subsection we calculate the derivative ^ 3 This is formally given 



dV\6 r ] 1 



dm 2 2 



<fx < (h(x) 2 > 



1 Those for the other theories can be derived similarly. 

2 Hereafter we omit the argument m . 

3 See the appendix in detail. 
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— / d 4 xd 4 y < <f){x)(j)(y) > 8{y - x) 



(4) 



d 4 x(j) c (x) 



d 4 xd 4 y 



<5 2 (-r[0 c ]) 



roc 
-oo 



5(j) c {x)8(j) c {y) t 

d 3 x" . We use this notation 



here "/ d 4 x n is an abbreviation of "Jo dr /f 
from now on. 

Since Eq.(|J) is a functional equation, it can not be solved directly. We, 
therefore, limit functional space and expand r[0 c ] in powers of derivatives [Bl, 
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d 4 x 



■7(^)-^W 



-X s (0 2 )(V0 c ) 2 + 



(5) 



where • • • denotes terms with higher derivative. Note that the coefficient 
functional of (do<f) c ) 2 differs from that of (V(/> c ) 2 due to absence of the 4- 
dimensional Euclidean symmetry. We then expand the both side of Eq.(|4]) 
with respect to derivative as Eq.(^) and match the coefficient functionals of 
each terms. In practice, we have to truncate the series in Eq.(^) at finite 
terms. In the present paper, we leave the three terms in Eq.(j5]). This is the 
next-to-leading approximation of the derivative expansion. We obtain the 
leading approximation, which corresponds to the evolution equation of the 
previous paper Reff28|, |29|j , by taking K = 1 and K s — 1 
From Eq.(||), l.h.s of Eq.(f|) becomes, 



dT 

dm 2 



d 4 x 



dV ldK Q 



2 + ^(V0 c ) 2 



From Eq.( 



M yx 



dm 2 2 dm 2 2 dm 2 

up to the second derivative terms,Q 

s ^ 2 r[0 c ] 

Scf) c (x)S(j) c (y) 
= 5{y-x)\v"{(j> c {y)) 



(6) 



4 Similar derivation in ref . |3^, |34[| has a problematic point: expiqx is handled as not a 
distribution but an ordinary function. 

5 Hereafter, we also use notations like V = V(4>) as values of functions at constant 
configuration, (p c (x) = 4>. 
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+K' ( ( f> c (y))d 2 oMy) + K Q ^ c {y)K) 

- {K K s , Oy <-> Vy}] . 

%-o;)(-F ^ -Z s V^+F") 

- {K «"> K s , dy Q «-> Vy}] 



mi 



(7) 



.4 



B„ 



A — 

B yx = 



yx ^yxi 

5(y-x)(-K d 2 yo -K s V 2 y + V") 

Ay X -Myx 

V"{<f> c (y))-V 

K (4> c (y)) - K 



(8) 
(9) 



Here, we divide M yx into two parts, A yx and B yx , which remains finite and 
vanishes at c (x) = <fi respectively. Details of the following calculation are ex- 
plained in the appendix. The right hand side of Eq. (|4]) can now be expanded 
around (j) c {x) = up to the second derivative as follows, 



d 4 xd 4 y( 



)- l 5{y-x) 



5cj) c (x)5(j) c (y)' 
d x(A 1 ) xx 

+ J dVyd^iA-^yByM' 1 )^ 

+ J d^xdSd^zd^v{A- l ) X yB yz {A- v ) zu B uv {A- l \ 



+ 



1 

px V 



1 

px V 



+ 1 -72 



--{<(9 O C ) 2 + <(V0 C ) 2 } 
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1 

px V 3 



+ 1 Ts 



2{K' (do<P c ) 2 + K' s (V<Pc) 2 }v' 



(10) 



px V 



-2 {K' p d (f) c + K' s pidi(j) c } {K p d (j)c + K s pidi(f) c } v' 
-\{v'f {F o (9o0 c ) 2 + F s (V0 c ) 2 } 



+ • • • + (terms which vanish at (f) c (x) = <fi like V x (something)), 

where v = K pl + K s p 2 + V", u' = K' p 2 + K' s p 2 + V'" (with p = 2imT(n = 

0, ±1, • ■ ■) and p 2 = Ef=i pf) and is an abbreviation of "T E£°=-oo / d 3 p/(2n) 3 J d A x" . 

We match the both side of Eq.(^j) using Eq.(|6|) and Eq.(JTD) and equate 
the coefficient functionals of (d<f) c ) respectively. After the matching, we put 
4> c (x) = 4> and get the following simultaneous partial differential equation, 



dV 
dm 2 
dKp 
dm 2 



dKs 
dm 2 



+ 



2 J P v 



p v A 

+k;i 4 + 2K n I - 



V 



V 



-2K K„ / pi- - -K a 



-X 

2 s 



n 2 



u 1 



p 

P 2 



+^ ^ + 2K. / — 



V 



—K S K S p^- - -K t 



n 2 



(11) 



(12) 



(13) 



where, U J " is an abbreviation of u TY,T=-oo I d p/{2ir) . Since we put 
4> c ( x ) — <P finally, we do not have contribution from the terms in Eq.(|T0|) 
which vanish at <p c (x) = <f). This is the evolution equation of the next-to- 
leading approximation of the derivative expansion. 
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Figure 1: Diagram of the external-momentum dependent self-energy at one- 
loop. 

2.2 Initial Condition 

We can calculate the initial condition M 2 ) using the perturbation theory 
within one-loop level thanks to the large auxiliary mass, M ~ T. The effec- 
tive potential, V(4>;M 2 ), is calculated within the one-loop approximation as 
follows, 

V(0; M 2 ) = ^M 2 + + 1 ^log(p 2 + p 2 + M 2 + ^0 2 ). (14) 

The loop correction to K (<f); M 2 ) and K s (4>; M 2 ) comes from the self- 
energy graph in Fig.[I] at one-loop level, which depends on external momen- 
tum, 



A 2 2 



n(g W^ 



1 



1 



2 J P p l + p 2 + m 2 + ^ (g - p ) 2 + (q - p) 2 + M 2 + ^ 



(15) 

The initial conditions of K Q ((p; M 2 ) and K s (<p;M 2 ) are given by the coeffi- 
cients of — q 2 . and — q 2 respectively. We, therefore, get 6 



K Q (4>\ M 2 ) 



dU(qlq 2 " 



1 - 



dqf 



o 



(16) 



<?i=0,<?o=0 

6 Strictly speaking, we first put qi = then put go = in taking the limit qi = and 
qo = 0. For the function Kq to be defined this order is essential, while for K s the order of 
taking the limit gives no difference. See ref. p5|. 
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Ks(b M 2 ) 



1 - 



1 + 



2 , p 
(in(g 2 , q 2 



(p 2 + p 2 + M 2 + 



(p 2 + p 2 + M 2 + ^f) 4 



<iq 2 







A0 2 



90=0, <ji=0 
1 



6 ^ ( p 2 + p2 + M 2 + A^)3' 
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We have calculated the evolution equation in §2^. and the initial condition 
2~2| . Since some of them, e.g.. the one-loop effective potential, has an ul- 
traviolet divergence, we have to renormalize it by the counter terms. Instead 
of considering this contribution, we simply assume that the renormalization 
effect is small and discard it. The contribution is, actually, small compar- 
ing with the finite temperature contribution around a critical temperature 
in most cases. We thus need to deal with only the temperature dependent 
pieces in the integrals of and the initial condition. To do so, we only 



perform the following replacement in Eq . (|Tl| 



with 



d 4 p 



1 



(27r)Hexp(p /T) 



We note that under the local potential approximation that K — K s — 
1, the evolution equation of the effective potential reproduce our previous 
equation. [BR BJI . 



3 Results and Summary 

We solved the simultaneous evolution equation Eq.([0]), Eq.([L2|) and Eq.(|T3|) 



with the initial condition Eq.(14), Eq.(lC) and Eq.fll7D numerically. We used 



an extended Crank-Nicholson method, which is explained in the appendix 
of Refp9|, to solve the partial differential equation. We could solve the 



equation at most temperatures above the critical one. Unfortunately, we can 
not, however, solve the equation at temperatures very close to the critical one 
due to numerical problems. We tried several improvements of the numerical 
method but failed. We, therefore, show only obtained results in the present 
paper. This equation may be solved at arbitrary temperature if excellent 
numerical methods to solve a partial differential equation are invented. Or 
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Figure 2: Effective potential around the critical temperature (A = 1). 

ability of computers progresses highly so that we can calculate with much 
higher precision. We use a mass unit, /i = 1 

We show the effective potential above the critical temperature in Fig|J 
for A = 1. We observe behavior of a second order phase transition up to 
this temperature. This effective potential seems to be cone-shaped. The 
critical temperature is estimated to be lower than that of the local potential 
approximation by 2%. 

We show K Q and K s in Fig.|3] for A = 1 . They are very similar in spite of 
the violation of the Lorenz invariance in this case while the initial condition 
is quite different. 

In summary, we derived an evolution equation of the effective action with 
respect to the mass squared in the ^-invariant scalar theory. We then ap- 
proximated the effective action by the derivative expansion. We showed that 
the previous evolution equation of the effective potential can be derived as 
the leading approximation, the local potential approximation. We next de- 
rived the evolution equation of the next-to-leading approximation which is 
the simultaneous partial differential equation. We finally solved the equation 
numerically. Though we could solve it at most temperatures above the crit- 
ical one, we could not do it under a certain temperature very close to the 
critical one unfortunately. However, this equation may be solved at arbitrary 
temperatures if excellent numerical methods to solve a partial differential 
equation are invented or ability of computers progresses highly. Anyway, we 
constructed the systematic method to improve the auxiliary mass method in 
principle. 

We are suported by JSPS. 
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Appendix 



In this appendix, we show the derivation of dT /dm 2 in detail. Some nota- 
tions are given in the chapter 2. The effective action T[cj) c ] for the Lagrangian 
Eq.(|l|) is defined as usual, 



(18) 



T[0 C ] EE W[J m ] - J rf 4 xJ m (x)0 c (x), 
W[J m ] EE log(Z[J m ]), 

Z[J m ) ee Jv[<t>] exp (| d 4 x£,E 



The equation(HI) is, then, derived as follows, 



dT[4> c 
dm 2 



Z[J n 



V\ 



,4 1 / ,/ n2\ i 9J m (x) 



Ix) > exp 



d 4 x£i 



,4 dJ m (x) , , 



/ d A x < Mx) 2 > 



d x 



(x) 2 - t j d'x{M)- l x . 



(19) 



Since Eq.(ffl) is a functional equation, it can not be solved directly. We, 



therefore, limit functional space and expand T[cj) c ] in powers of derivatives^ 
(seeEq.©) 



n 



d x 



-V(€) - \K {<p 2 c ){d Q <p c ) 2 - ^ s (0D(V0 c ) 2 + ■ ■ 



where "• • ■" denotes terms with higher derivative which are omitted here and 
equating the coefficient functionals of (d<p c ) in the both side of Eq.(f|). 

The left hand side of Eq.([Tj]) is calculated to be Eq.(|B|) by simply differ- 
entiating T with respect to m 2 , 



dT 

dm 2 



d 4 x 



dV ldK ldK 2 



dm 2 2 dm 2 



2 dm 2 
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It is very complicated to calculate r.h.s. of Eq.([T^). First we calculate 
M yx . The derivative of the effective action T[(f> c ] with respect to c (x) is, 

Sir ~ <20) 

+^(0 c (x))(«9o0 c (x)) 2 + Ko(0 c (x))(d o 2 c (x)) 
+ {K ^K s ,d ^ V}. 

Here, we assume that we can make a partial integral freely without a surface 
term. 

We define the operator M xy in the following sense. An operator O xy de- 
fined through functional derivative with respect to (p(x), say 5F(<ft(y))/5<j)(x), 
acts on any appropriate test function, say T(y), 

OT{x) = [d 4 y{-J—F(^y)))T(y). 



5<f)(x) 

In particular, if F(<f>(y)) contains the derivative of (f)(y), say F(4>(y)) = 
G(<t>(y))d<f>(y), 

OT(x) = J d 4 y S.jJ^G((j)(y))d(j)(y)^T(y). 

= [ d A y8(y - x) [G'(0(y))^(y)T(y) - d {G(0(y))T(y)}l .(21) 



In order to obtain Eq.(^) and to determine the inverse of M yx around 
we divide the Eq.(|7|) into two pieces (Eq.(||)), 



M-yX Ay X By X , 



A yx = 6(y-x)(-K d 2 -K s V 2 + V ), 
By* = -5(y-x)[v"(<f> c (y)) 

(~<(<^)) W> c (y)) 2 + K' (<f> c (y))d o( f> c (y)d yo 



.2 

+ K' (<p c (y))d 2 oMy) + Ko(Mv)K} 
- {^o <-> k„ d yo <-> vj] . 



n 



The inverse of M yx is, then, expanded as, 

n=oo 

= {(A-^E^DW (22) 

n=0 

Here, the multiplication of the "matrix" A and B is taken in the following 

sense, 

{AB) xy = J d 4 zA xz B zy . 
The inverse of A xy is calculated easily, 

(A- 1 ) xy = [ -esxp{ip{x-y)}, (23) 

Jp Up 

where v v = K pl + K s p 2 + V . 

We need terms with (d(f> c ) n (n = 0,2) to evaluate r.h.s of Eq. (|19"D . Such 
terms are contained only in the first three terms in Eq . (|2~2"D , 

(M-% y = (A~% y 

+ (A- 1 BA- 1 ) xy (24) 
+(A- 1 BA- 1 BA- 1 ) xy . 

The relevant terms in Eq.(|l^) are, therefore, the following, 

Jd'xiM- 1 )^ = Jd^A- 1 )^ (25) 
+ J d i x(A~ 1 BA~ 1 ) xx (26) 



+ J d A x(A- l BA- 1 BA- 1 ) xx . (27) 



We next calculate Eq.(p5|)- (|27|) up to terms with the second derivative of 



x). From Eq.(|2"3"|) the first term (|2li|) is easy to calculate 



d A x(A~ l ) xx = I d 4 x f -. (28) 
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The second term Eq.(p6|) is 
J d 4 x{A- l BA~ l ) xx 
= J d 4 x J d 4 y J d 4 z J — exp {ip(x — y)} 

5{y-z) 

+ K'MyMMy) + Ko(My)K} + » <- k„ d yo <- vj 

/ — exp {iq(z — x)} 

J a Vn 



v"(Uv)) + [-K(MyW My)) 2 + <(0 c (y))d o c (y)a, 



hi i 



(integrating over z.) 

[ d A x I d 4 y I I exp {ip(x — y)} exp {iq(y — x)} 

J J Jp Jq Vp Vq 

■v"(Mv)) + {^K{Uy)){d^c{y)f + <(</> c (y)W c (</)% 

+ K' Q {<t> c (y))dl<t> c {y) - MMy))q 2 o} + {Ko <- K s , d <- V, q <- q} 
(integrating over x and q then replacing y by x.) 
1 



d"x 



7 2 
P 



-V (cf> c (x)) + j-# (0 c (x))(<9 o 0(x)) 2 + K Q (<p c (x))d <P(x)i Po 

+ K'^x^dl^x) - k (<f) c (x))pl} + {K <- K„ d <- V,p <-> p} 
(partially integrating the fourth term.) 



d 4 x I — 

p K 



-V"{Ux)) ~ { -A^(0 c (x))(9 o c (x)) 2 +K ((f )c (xM 



- {K K s , d <-» V,p <-> p}] • 



rf 4 x / - 



~<(d O c ) 2 + V s '(V0 c ) 2 



+ (ierms which vanish at c (x) = Zz&e V" x (something)) . (29) 

Since we set c (x) = after the matching of Eq.(|J), we have no contribution 
to the evolution equation from the terms which vanish at (j) c (x) = 0. We 
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thus leave only terms which remain finite at (j) c {x) = <fi from now on. The 
third term fl27D is very complicated to evaluate, 

J d^xiA^BA^BA- 1 )^ 
= J d 4 x J d 4 y J d^z J d A u J d 4 v J — exp {ip(x — y)} 

S(y - z) [-V"(&(y)) + S^K'^UyWoMy)) 2 + K' (<p c (y))d <p c (y)dyo 
+ K' (<f> c (y))d 2 <f> c (y) + K (My)K} + {K <- K s , d yo <- VJ 
/ — exp {iq(z — u)} 

JqVq 



5{u — v) 



\ ' (o,.( „)) + { ^K;(<p c (u))(d <p c (u)) 2 + K' {<j> c {u))d Q <j> c {uy), 



<<<> 



+ K' (<p c {u))dl<p c (u) + k (<p c (u))d 2 U0 } + {K <- K„ d uo <- VJ 
exp {ir(v — x)} 



r V r 



(integrating over x, z, v, r and replacing u with x) 
f d 4 x f d A y f f -^—exp{i(p-q)(x-y)} 

J J JpJq V* U q 

-V"(<f> c {x)) + ^-K'^ c {x)){dMx)) 2 + K' {(f) c {x))d 0( f> c {x)ipo 

+ K' (<j) c (x))d 2 o( j) c (x) - K (cf> c (x))p 2 } + {K <-> K a , d <-> V,p <-> p} 
1 



+ 



■v"{My)) + [^K{<p c {y))(dMy)) 2 + K' (<p c (y))doMy)m 
KiMvMMy) - MMy))ql} + {k <- k„ d <- v, qo <- q} 



(partially integrating the term with d^cf).) 
f d 4 x f d A y f f \—exp{i(p-q)(x-y)} 

J J JpJq Vt V n 



ipJq ,s p u q 
1 



-V {<j> c {x)) + {--K (cj> c (x))(do(i> c (x)) ' + K (<p c (x))dMx)iq 
- K (4> c (x))p 2 } + {K *-* K s , d <-+ V,p *-* p, g <-> q} 

-v"(My)) + (MymdoMy)) 2 + KiMyVdoMv)™ 



14 



- K ((f) c (y))ql} + {K <-> K 8 , d <->V,p ^ p, go <-> q} 
By the following variable exchange, 



P 
5 



p + g 
P 



we obtain, 

= f d A x I d A y I I — - exp {iq(x - y)} 

J J Jp Jq Vp+q Vp 



X 



-V (0 c (x)) + |--K o (0 c ( a :))( ( 9 o c ( a :)) i + K o (0 c ( a :)) ( 9o0 c (x^ o 
- K ((f) c (x))(p + q)lj + {K <-> K s ,d <-»• V,p <-> P, go <-> q} 
-^"(0c(y)) + |-i<(0 c (i/))(9o0 c (|/)) 2 + <(0 c (y))9 o c (y)i(p + g) 



- K (<f) c (y))pl} + {K <-> K s , d <-> V,p <-> P, go <-> q} 



P J 9 



{^o(4p go + 2g 2 ) + i^(4pq + 2q 2 } 12 {x Pogo + ^pq} 
< 1 - H s - 



+ 



/A, 



-^"(0c(x)) 



•2 i ^o( < / , c(a;))(5o0 c (a;)) 2 + X o (0 c (z))d o c (a;)ipo 



x 



- K ((f) c (x))(p + q)lj + {fT 4-* K s ,9 <-> V,p <-+ P,go <-> 
-t/"(0 c (y)) + Wc(l/)) 2 + <(0 c (y))9 o c (y)i(p + g) 



- K o (0 c (y))p 2 } + {K <- tf„ 9 ^ V , Po <-»• p, go <-»• q} 

where we take the terms up to g 2 since we will replace g to the derivative and 
need terms with second derivatives. Similarly, by the following exchange of 
the variables, 



p — ► p 

q — ► p + g ' 
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we have, 



d 4 x ( d 4 y I I — — \ exp {iq(x - y)} 

J Jp Jq V p +q V p 



-V\<t> c {x)) + |-i<(0 c (x))(9 o c (x)) 2 + K' (<f> c (x))d ci> c (x)ip 

K ((f)c(x))(p + q)l] + {K «-> K s , d «-> V,p «-> P, go <-> q} 

-y"(0 c (y)) + (MvWoMy)) 2 + K' {<p c (y))do<Pc{y)i{p + q)o 

Ko{<Pc{y))p 2 } + {K «-> K s , d ^V,p ^ p, g «-> q} 
exp {ig(x-y)} 



d x d y 



pJq 



(31) 



^ {A (2p g + g 2 ) + ^s(2pq + q 2 )} 4 {A> g + A s pq} 



+ 



1 



-y"(0 c (x)) + |-|A o '(0 c (x))(9 o c (x)) 2 + K' Q {<l) c {x))do<t> c {x)ipo 

- K ((f)c(x))(p + g)o} + {A <-> A s ,<9 «-> V,p «-> P,go <-> q} 
-y"(0 c (|/)) + {-i<(</> c (y))(a o c (y)) 2 + KiUv^Uv)^ + g) 

- A o (0 c (y))Po} + {^o <-> a s , <9 «-> v,p <-> p, go <-> q} • 

By comparing Eq.(|30f) with Eq. fl3lD , we have the following "identity", 
exp {iq(x — y)} 



d x Id y 



pJq 



(32) 



{K (2p qo + ql) + K s (2pq + q 2 )} 8 {K p qo + A s pq}' 



-V"(<p c (x)) + { --K' \<p c (x))(do<P c (x)) 2 + K' (<p c (x))d Q <P c (x)tpo 



- K (4> c (x))(p + g) 2 ,} + {A <-> A s ,<9 <-> V,p <-> P, go <-> q} 

-y"(0 c (i/)) + {-^<(0 c (2/))(a o c (y)) 2 + KttMfaMv)* 

- Ko(4> c (y))pl\ + i K o *-> #o <-> V,p <-> p, go <-> q} 
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Indeed, this identity holds exactly for the coefficient functions of (V0 C ) 2 , 
though it is not the case for that of (<9o0 c ) 2 . Anyway, we use this identity 
and rewrite (or define) the third term Q27D in the following, 

J d^xiA^BA^BA- 1 )^ 

^/A// eXPfa( 3 X " y)} (33) 

J Jp Jq U p 

^ 1 {K {2p qo + ql) + ^ s (2pq + q 2 )} 

2 u p 

-V'\4> c {x)) + |-i<(0 c (x))(9 o c (x)) 2 + K' (<p c (x))doMx)Wo 

- K {(f) c (x))(p + q)lj + {K Q ^ K s ,d «-> V,p «-> p, q «-> q} 

x -v"(4> c { y )) + {-ii^(0 e (y))(^c(y)) 2 + KiMvVdoMvW + ?)o 

- ^o(0 c (y))Po} + W <9o <-> V, J9 p, g <-> q} • 

Note that there is no contribution from Eq. ( |3"3"| ) to the evolution equation for 
the effective potential V since the constant terms always include K and V" 
which are zero at <p c (x) = (f> by definition. 

From now on we evaluate the coefficient functions of (d<f) c (x)) 2 . First we 
calculate the first term of Eq.(|33]). There are three contributions from this 
term, 



d A x I ' d A y I I 

J Jp Jq 



exp {iq(x-y)} 



[K' ((f) c (x))d (f) c (x)ip + K' s ((f) c (x))V(f) c (x)ip} 

x [K' ((j) c (y))d Q (j) c {y)ip + K' s {<t> c {y))V<j) c {y)iv) 
(integrating over q and y.) 

= fd 4 x / 4 {Kido^Po + ^l(V0 c )p} 2 , (34) 



d 4 x I ' d 4 y f f 

J Jp Jq 



exp {iq(x — y)} 
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[K' ((j) c (x))do(j)c(x)iqo + K' s ((f) c (x))V(f) c (x)iq} x 

{-MMvM - K s (<p c {y)) V 2 - V"{<j> c {y))} 
+ 2 {K' (4> c (x))d (f) c (x)ip + K(0 c ( a; )) v 0c(a;)«p} x 

{-K (My))poqo - ^ s (0 c (z/))pq} 

'{K' {(f> c {x))d <f> c {x)d X0 + K' s ((j) c (x))V(j) c (x)V x } x 

{-K o (0 c (y)K - K s (<p c (y))p 2 - v"(My)} 

+ 2 {Ko(0 c (x))<9o0 c (x)ipo + ^(0c(x))V0 c (x)ip} x 
{^o(0c(2/))iPo9a: o + K s ((f) c (y))ipV x } 
x exp — y)} 

(partially integrating over x and then integrating over q and y .) 



d x 



1 r 



{<(a o c ) 2 + <(V0 C ) 2 } {<p 2 + <p 2 + V'"} (35) 



+2{<(9 o 0>o + <(V0 c )p} i 



and 



d 4 x / d 4 ?/ 



exp {ig(x - y)} 



p J q 



{-^o(0 c (x))g o 2 -^ s (0 c (x))q 2 }x 

{-iro(0 c (2/))Po " K s {<f> c (y))p 2 - V"{<f> c (y))} 



d x 



1 r^/ 



{^(9 O C ) 2 + <(V0 C ) 2 } {<p 2 + <p 2 + V'"} . (36) 



As a whole, we have the following contribution from the first term of Eq. 

1 



d"x — n 



2{K' (d 0( p c ) 2 + K' s (V ( p c ) 2 }v' p 



(37) 



+ {K' d (f) c p + <V0 cP }' 
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Similarly, we have the following contribution from the second term of Eq. (|33|) , 

-2 [K' (d <P c )po + ^l(V0 c )p} [K (do<Pc)po + K s (V<p c )p} v' p ] . (38) 

Combining Eq. fl28|) , Eq. (|29"D , Eq.(|37D and Eq. (|38|) , we obtain the evolution 
equation, Eq. flToD . 
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